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Abstract 

We apply the linear (^-expansion (LDE), originally developed as a nonpertur- 
bative, analytical approximation scheme in quantum field theory, to problems 
involving noninteracting electrons in disordered solids. The initial idea that 
the LDE method might be applicable to disorder is suggested by the resem- 
blance of the supersymmetric field theory formalism for quantities such as 
the disorder-averaged density of states and conductance to the path integral 
expressions for the n-point functions of Xcj) 4 field theory, where the LDE has 
proved a successful method of approximation. The field theories relevant for 
disorder have several unusual features which have not been considered before, 
however, such as anticommuting fields with Faddeev-Popov (FP) rather than 
Dirac-type kinetic energy terms, imaginary couplings and Minkowskian field 
coordinate metric. Nevertheless we show that the LDE method can be suc- 
cessfully generalized to such field systems. As a first test of the method and 
also to give some understanding of its origins, we calculate to third order in 
the LDE the ground state energy of a supersymmetric anharmonic oscillator 
with FP kinetic term and real anharmonic coupling strength of arbitrary mag- 
nitude. Strong evidence for the convergence of the LDE is obtained. We then 
calculate to second order in the LDE the disorder-averaged density of states 
of a one dimensional system and find even at first order more accurate results 
than the commonly used self-consistent Born approximation. In the final 
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part we outline one possible way in which the LDE method might be applied 
to the conductance, using as supporting example a zero-dimensional model 
with Minkowskian field coordinate metric. Further directions for research are 
discussed in the conclusion. 

PACS numbers: 02.60.Gf, 11.15.Tk, 71.23.-k, 72.80.Ng 
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I. INTRODUCTION 



Interacting quantum field systems rarely allow exact, closed-form expressions for their 
observable quantities. In nearly every case some method of approximate solution must 
be employed. The most common approximation involves developing a series expansion 
of the quantity of interest in increasing powers of the coupling constant multiplying the 
interaction term in the field Lagrangian. However, such a series is not expected to give a 
good approximation when the coupling constant in the appropriate dimensionless units is not 
small. Even when the coupling is small, there are situations where the series approximation 
does not work, such as when the quantity of interest is nonanalytic at zero coupling. Given 
the relevance of strongly interacting quantum field systems for describing processes occurring 
in nature, a worthwhile and stimulating challenge has been to find alternative, analytical 
approximation schemes which must necessarily be nonperturbative in the physical coupling 
constants of the field systems. 

The linear 5-expansion (LDE) is one such scheme which has been successfully applied to 
problems in, for example, A0 4 theory0 and quantum chromodynamics.ii In outline, the LDE 
method replaces the action S of the field theory with a modified action S$ which interpolates 
linearly between a soluble action So depending on a variational parameter Q and the original 
S, i.e., 



The Green function of interest G is evaluated as a power series in the artificial parameter 5 
up to the desired order iV and then 5 is set equal to 1. Unlike the exact Green function G, 
the truncated Green function Gn will depend on Q. The latter is fixed at each order N by 
applying the principle of minimal sensitivity (PMS) J 



The PMS condition is crucial, in that it provides the nonperturbative dependence on the 
physical coupling parameter and in several cases has been shown to ensure convergence of the 




(1) 
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sequence of approximants Gn(&n)- There is considerable arbitrariness in the choice of Sq. 
Some choices will result in more rapid convergence of the sequence or simpler calculations 
to perform, whereas other choices may not work at all. At the very least, Sq should have 
some resemblance to the original action S. Proofs of the convergence of the LDE method 
have been obtained in several different ways for certain Green functions of the quantum 
anharmonic oscillator.H These proofs have also given us insights into how the method 
works. For higher dimensional field theories no convergence proofs have been constructed 
so far and we must rely on comparison with numerical methods. 

In this paper we apply the LDE method to a different problem: the quantum dynamics 
of a single electron moving in a random background potential. The idea of using the LDE 
method comes from the so-called supersymmetric field formulation of the problem,§0 in 
which quantities characterizing the quantum dynamics, such as the disorder-averaged den- 
sity of states and conductance, are represented as path integrals over both commuting and 
anticommuting field coordinates of some interacting field system. For a random potential 
which is Gaussian delta-function correlated the corresponding supersymmetric field actions 
resemble the usual A0 4 action, hence suggesting the possibility of applying the LDE method 
to the former systems as well. Of course, there are already several well-established approxi- 
mation methods for studying the quantum dynamics of an electron in a random potential. 
Nevertheless, we thought it would be interesting to see whether a method originally de- 
veloped for the study of relativistic quantum field systems could be successfully adapted to 
random electron systems. First indications are that the LDE method in fact performs rather 
well in comparison with established approximation methods. 

The supersymmetric field actions have several unusual features which have not been 
considered before in LDE investigations. The kinetic term for the anticommuting fields is of 
Faddeev-Popov (FP) rather than the usual Dirac form (i.e., d t x*d t x instead of x*dtx)- Also, 
the interaction terms appearing in the actions are imaginary and, for the action associated 
with the disorder-averaged conductance, we have a Minkowskian field metric [i.e., r\ mn 
with r] mn = (—l) m S mn ; m, n = 1, 2]. In the following sections we show how the LDE method 
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generalizes to such systems. 

In Sec. II we evaluate to third order in 5 the ground state energy of a supersymmetric an- 
harmonic oscillator with FP kinetic term and real, positive anharmonic coupling parameter. 
The action for this system is in fact the same as that associated with the one-dimensional 
density of states for delta-function-correlated disorder, the only difference being that the 
coupling is real instead of imaginary. Although the ground state energy is of little relevance 
for the random electron dynamics, it is the simplest quantity on which to test the LDE 
method. Furthermore, the calculation gives some understanding of the method's origins by 
showing that to first order in 5 it is equivalent to choosing a Gaussian wavefunction with 
variational parameters which are fixed by requiring that the energy expectation value be 
a minimum. As we shall see, the LDE method reveals a rather remarkable dependence of 
the ground state energy on the frequency and anharmonic coupling parameters, making the 
supersymmetric anharmonic oscillator an interesting system in its own right. 

The LDE method is applied to the disorder-averaged energy density of states in Sec. III. 
We first consider a "zero-dimensional" density of states model (i.e., path integral replaced 
by ordinary integral), for which we can easily go to high order in S and obtain support 
for the convergence of the expansion. We then evaluate to second order in 5 the one- 
dimensional density of states for delta-function-correlated disorder and compare with the 
exact expression, as well as with the self-consistent Born approximation (SCBA). While 
the use of the LDE was motivated by the supersymmetric field formulation, it is in fact not 
necessary to calculate the density of states within this formulation. Indeed, we show that the 
same series approximation is obtained by expanding in 5 the averaged one-electron Green 
function expressed directly in terms of a ^-modified single-electron Hamiltonian. Note that 
the LDE method is not restricted to low dimensions, nor to just delta-function-correlated 
random potentials. 

Because the supersymmetric action associated with the disorder-averaged conductance 
has Minkowskian field coordinate metric, the action must be 5-modified in a way which is 
quite different from that of the density of states action. In Sec. IV, we describe one possible 
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^-modification. In contrast with the density of states, the supersymmetric field formulation 
is essential for the application of the LDE method to the conductance: there is no modified 
physical Hamiltonian which yields the same series expansion in 5. As a first test, we apply 
the method to a zero-dimensional conductance model. 

In the conclusion we suggest further directions for research. This includes the application 
of the LDE method to the supersymmetric non-linear sigma model in order to determine 
the conductivity critical exponent. 



II. THE SUPERSYMMETRIC ANHARMONIC OSCILLATOR 

In this section we use the LDE method to determine the ground state energy of a quan- 
tum, supersymmetric anharmonic oscillator. The classical Lagrangian is 

L = i$t$ _ ^$t$ _ A($t$)2 ; ( 3 ) 
2 2 



where 



$t = ( x * s*) (4) 



is a supervector coordinate with anticommuting component \ and commuting component 
s. In all our calculations involving anticommuting variables we follow the rules and conven- 
tions of Efetov,0 the only difference occurring in the definition of complex conjugation for 
anticommuting variables. Our definition is 

(X1X2T = Xlxl, (**)* = X- (5) 

The above Lagrangian is motivated by the supersymmetric formulation of the disorder- 
averaged density of states (see Sec. III). In one dimension and for delta-function-correlated 
random potential, the supersymmetric Lagrangian associated with the density of states 
differs from Eq. ([3|) only in having an imaginary instead of real coupling (as well as an 
unimportant overall sign and constant coefficients). Bohr and Efetovffl also investigated 



the above supersymmetric Lagrangian and found a closed form expression for an eigenstate 
with eigenvalue zero which they assumed to be the ground state. As we shall see, however, 
this assumption is correct only for u 2 larger than a certain negative value. 

In order to quantize the oscillator we must first write down the Hamiltonian. In terms 
of real coordinates S\, s 2 , Xi> an d X2, where s = Si + is 2 and x = Xi + the Hamiltonian 
is 

H = SiPi + XiTd ~ L 

= \{Pi + Pi) + 7^ 2 ( s i + S D - "Ti*2 + iuj 2 XiX2 + K 2i XiX2 + sl + si) 2 , (6) 

where the momenta are defined as follows: 

dL dL 
n = W = °* ^=g~ = ^Xr (7) 

The oscillator can now be straightforwardly quantized using the correspondence principle, 
i.e., by associating with the position coordinates and conjugate momenta operators acting 
on some state space and satisfying canonical (anti) commutation relations: 

Sj) = -iSij and {hi, Xj) = -i<% , (8) 

with all other (anti) commutators vanishing and where we have set % = 1. Note from Eq. 



(0) that 7Tj is antihermitian, i.e., ftj = — 



7T; 



It will be convenient to work with two different state space bases. In one, the states are 
wavefunctions depending on the position coordinates, , s 2 , Xi> X2), with scalar product 
defined as 

= ~i J dxidx2ds 1 ds 2 ip*(si, s 2 , xi, Xa)^(si, s 2 , Xi, X2), (9) 

and the position and momentum operators represented as follows: 

d d 
Si ^ s u pi <r-> -i—, Xi «-> Xi, «-> -i^—. (10) 
OSi oxi 

The energy eigenstates of the free Hamiltonian (A = and u 2 > 0) form a second useful 
basis. This basis is most easily constructed using creation/annihilation operators, defined 
in terms of the position and momentum operators as follows: 



2u 

Pi = _i ^( a _ flt ) 



S-2 



' (6 + 6 f ) 



2u 



l>2 -'Xl^ip-ti) 



Xi 



7T1 = -. 2 



-^(c + c 1 " + d + d f ) 



i , 



X2 = ~ l 7T~r^[ c ~ c ' — d + (tf) 



,t 

'2^ 



7T 2 = -^( C - C t + d_^), (11) 

where we have omitted the hats on the operators. From these definitions and the canonical 
(anti) commutation relations @, we have 



a, a 



1, [6,6+1 = 1, {c,ct} = l, {d,rft} = -l, (12) 



with all other (anti) commutators vanishing. Using Eqs. (O) and fll2"D to express the non- 
interacting part of the Hamiltonian in terms of the creation/annihilation operators, we find 
that 

H = u{a ] a + b ] b + c ] c- d ] d). (13) 

The energy eigenstates are obtained in the usual way. We first introduce a state |0), satisfying 

a\0) = 6|0) = c|0) = d\0) = (14) 

and (0 1 0) = 1. ^From Eqs. fll3D and (0) we see that |0) is an energy eigenstate with 
eigenvalue 0. All other eigenstates are obtained by acting on |0) with the creation operators 
a), ftt, c\ and d) . The eigenstate 

K , iiH ffifflW m (15) 

has eigenvalue 



Emnjk = {m + n + j + k)u > 0, (16) 

for integers m, n, j, and k > 0, so that |0) is in fact the ground state of the free oscil- 
lator. Note that, since c' and d) anticommute, j and k are either zero or one. Note also 
that, because of the minus sign entering in the anticommutation relation (|l|) for d and d\ 
eigenstates with k = 1 have negative norm. The occurrence of negative norm states in the 
state space can be traced back to having an anticommuting coordinate kinetic term with 
two time derivatives, x*x, in Eq. (0). The FP ghost fields, which occur in the quantization 
of Yang- Mills fields, have similar properties .111 Negative norm states can be avoided by using 
instead a Dirac kinetic term, x*X- However, since the FP kinetic term is the relevant one 
for disorder and the Lagrangian (^) is not meant to describe a physical oscillator, we shall 
live with the negative norm states. 

Using Eqs. (|l4|), and (|To| ) we find that the ground state for the free Hamiltonian 



has the following form in the position coordinate basis: 



^oOi,S2,Xi,X2) = ^^ ex P P ^ t$ ) = 7^ exp ~\ UJ ( s l + s l + 2i XiX2) • (17) 

This wavefunction is just the supersymmetric generalization of the Gaussian function. 

Let us now return to the problem of determining the ground state energy of the su- 
persymmetric oscillator Hamiltonian for A > 0. In the case of the ordinary, commuting 
coordinate anharmonic oscillator it is well known that the minimum energy expectation 
value for a Gaussian wavefunction with variable frequency and coordinate shift parameters 
gives a good approximation to the ground energy.O A natural choice for the supersymmetric 
oscillator trial wavefunction is then 

1 /fi^A 1 / 4 



ip{si,s 2 ,Xi,X2) = -7=1 — —) exp 

V27T V V* J 



~fii(si - s ) 2 - ^2-4 - ivXiXs 



(18) 



where sq is the shift parameter and Q%, and v are the frequency parameters. A nonzero 
shift parameter is necessary in order to approximate the ground energy for u 2 < and 
large, since in this case the commuting coordinate part of the potential in (|B|) has the 



\0J 



"Mexican hat" form and the ground-state wavefunction peaks at the rim of the hat, instead 



of at the center. The expectation value (ip\H\ip) is most easily worked out in the position co- 
ordinate representation using Eqs. (^j), flip]), (IS), and (^). The minimum expectation value 



with respect to variations in the parameters sq, Q\, Q2, and v then gives an approximation 
to the ground state energy. 

The above method provides only a one-off approximation, however. A considerable 
improvement would be a scheme in which the expectation value (ij)\H\ij)), for given by 
Eq. fll8f) , appears as the first-order term in a perturbation series. As shown by Stevenson 
for the ordinary anharmonic oscillator (see Sect. V of Ref. [13]), such a scheme can in fact be 
realized. Generalizing to the supersymmetric oscillator, the scheme involves first modifying 
Hamiltonian (||) as follows: 

H s = 2(^1 + Pi) + 2~ fi i s i + \ Q 2 s l ~ i7ri7r2 + iy2 XiX2 + \ u2s l + As o 



+5 ^ 2 {s\ + s ) 2 + \u 2 s\ + iu 2 XiX2 + A (2^X1X2 + (si + s ) 2 + sf 

L 2 2 \ 4 1 r>2 2 - 1 - n2 2 ■ 2 



(19) 



and then solving for the lowest energy eigenvalue of this modified Hamiltonian using the usual 
Rayleigh-Schrodinger (RS) perturbation procedure with 5 serving as expansion parameter. 
Note that H$ interpolates linearly between a harmonic oscillator Hamiltonian for 5 = and 
the original anharmonic Hamiltonian (||) with coordinate redefinition si — ► si + So for 5 = 1. 
Hence the name: "Linear Delta Expansion" . To first order in 5, the RS series approximation 
to the ground energy is (with 6 set equal to 1): 

E w 4 0) + (0|#int|0>, (20) 

where E^ and |0) are the ground energy eigenvalue and eigenstate, respectively, of the 
harmonic oscillator Hamiltonian H$ =0 and H int = H$=i — Hs=o- The most straightforward 
way to work out Eq. ( PD| ) is to first express Hs in terms of creation/annihilation operators 
using relations ( |TTD with the appropriate frequency changes. After some calculation we find 
that Eq. (120) indeed coincides with the expectation value (ip\H\ip). However, we now have a 



systematic procedure which allows us to go beyond the Gaussian variational approximation: 
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Eq for the Hamiltonian H$ is evaluated up to the desired order N in 5 using the RS per- 
turbation method and the frequency and shift parameters are then fixed by minimizing the 
order N approximation to Eq. The calculations proceed in much the same way as for the 
ordinary anharmonic oscillator, the main difference arising from the negative norm eigen- 
states in the state sums which appear at second order and higher. The negative norm states 
are best dealt with by making explicit the eigenstate normalization factors in the RS series 
expansion formula. 

In Fig. [I] we plot the results of the order S and S 3 ground state energy calculations. [To 
order 5 2 , there were no parameter values for which Eq was stationary. This is a common 
occurrence in LDE calculations, where stationary points may exist for odd (even) orders 
only] The analogous results for the ordinary anharmonic oscillator are given in Ref. [13|. For 
the range u /A 2 / 3 > -3.69, the ground energy is exactly zero to order 5 3 . The frequency 
parameters satisfy £Xi = Q2 = is, while the shift parameter so = 0. For cu 2 /\ 2 ^ 3 < —3.69, Eq 
is stationary for a choice of parameters with the same pattern as above, again giving Eq = 0. 
However, there is another choice of parameters having the pattern 7^ Q 2 — v and Sq 7^ 
for which Eq is also stationary but less than zero, hence providing a closer approximation 
to the ground state energy. 

An independent check of the results for uj 2 < and \u 2 \ large is obtained by expressing 
the Hamiltonian @ in polar coordinates and expanding the commuting coordinate potential 
to quadratic order in the radial coordinate difference r — r , where r is the location of the 
potential minimum. To leading order, the ground state energy of the resulting harmonic 
oscillator Hamiltonian is 



a; 4 



(21) 



uj 2 — > —00. 



16A 

We have verified that the order 5 and 5 3 approximations indeed tend to fl2~T|) as 
While not constituting a proof, this result taken together with the closeness of the order S 
and S 3 approximations strongly suggests that the (^-expansion converges. 

The existence of a region where the ground energy is exactly zero is reminiscent of 
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quantum mechanical systems with time translation supersymmetry [i.e., the square of the 
supersymmetry operator equals the Hamiltonian (for a review, see e.g., Ref. |T^)]. For such 
systems, invariance of the ground state under supersymmetry transformations implies the 
vanishing of the ground state energy. On the other hand, if the ground state is noninvari- 
ant, time translation supersymmetry is spontaneously broken and we have a nonvanishing 
(positive) ground state energy. It is tempting, therefore, to speculate that time translation 
supersymmetry is broken for uo 2 below — 3.69A 2//3 and unbroken above this value. However, 
the Hamiltonian (f|) is known only to have rotational supersymmetry (i.e., the square of the 
supersymmetry operator equals the angular momentum operator corresponding to rotations 
in the s\— S2 coordinate plane) and no conclusions can be drawn concerning the ground en- 
ergy from the invariance properties of the ground state under rotational supersymmetry. It 
would be of interest to try to understand the reasons for the rather remarkable dependence 
of the ground state energy on uj 2 shown in Fig. [I]. 

III. THE DENSITY OF STATES 

Consider an electron in a d- dimensional random potential, described by the Hamiltonian 

H= m +V (r)-E, i = l,...,d, (22) 
2m 

where E is the Fermi energy and the potential V is Gaussian-distributed with delta-function 
correlation: 



V{r) = 0, V(r)V(r') = \5(r-r'). (23) 

The overline denotes disorder averaging and the parameter A characterizes the strength of 
the disorder. The supersymmetric formulation of the averaged energy density of states per 
unit volume for this system i 

pjE) = Im G(+)(r,r; E) = -Re / D&D® s(r)s*(r) exp(zS), (24) 

7T TT J 
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where is the retarded one-electron Green function. The supervector coordinate $ is 
defined in Eq. (Eh and the supersymmetric action S is defined as follows: 



S = - [dr —di&di® -(E + ie)&® - — (V$V 



(25) 



Given the resemblance of this action to the ordinary A0 4 action and the effectiveness of the 
LDE approximation method for studying various quantum properties of A0 4 theory,0 it is 
natural to apply the LDE method to the density of states of system ( p2]) as well. Of course, 
as we have already mentioned, the supersymmetric action has several significant additional 
features. In the preceding section we showed how the LDE method could successfully ac- 
commodate some of these features in ground state energy calculations. In the following, 
we further extend the LDE method, using it to approximate the expression (|24l) which is 
essentially the two-point function of a supersymmetric A0 4 system with imaginary coupling. 



The first step in the LDE procedure is to "5-modify" the action (|25|). Recall that, for 



the parameters satisfying Vt\ = Vt 2 = v and Sq = 0, the 5- modified Hamiltonian fllTf ) gave 



the energy Eq = 0. Furthermore, Bohr and EfetovS (see also Sect. 6 of Ref. [T(]) showed, by 
expressing the d = 1 density of states in terms of the eigenstates of Hamiltonian (^h, that 
only the eigenstate with eigenvalue zero was relevant. Noting the form of the 5- modified 
Hamiltonian (|I^) for Qi = Q 2 — v — ^ and So = 0, we are led to consider the following 
5-modified action: 



r i fi 2 

S s = - dvl — t $t0 $ _ (q + ^ e )$t$ + s 
J 2m 



(26) 



The original action ( p5|) is replaced by the modified action (p6|) in Eq. (|2^) and the series 
expansion in 5 obtained. The final step is the application of the PMS condition (|2]) in order 
to fix the frequency parameter Q. Note that the PMS condition is applied before taking the 
real part in expression (p4|). 

Again it is not necessary to implement the LDE method within the supersymmetric for- 
mulation of the density of states. Comparing the actions S and S$ and examining also the 
correspondence between the action S and the physical Hamiltonian ( f22|) with correlation re- 
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lation (^), we can immediately write down the modified physical Hamiltonian corresponding 
to the action S5: 



PiPi 
2m 



tt + 5 (n-E) + S 1/2 V(r). 



(27) 



Since the averaged potential satisfies V(r) = 0, only integer powers in S appear in the 
expansion. If we then replace the physical Hamiltonian by this modified Hamiltonian in the 
definition of the retarded one-electron Green function: 



G ( s + \r,r') = (r\(-H s + ze)- i \r'} 



(28) 



we recover the same series expansion in 5 as for the above supersymmetric formulation. 
Although the supersymmetric formulation is not required in order to apply the LDE ap- 
proximation method to the density of states, the formulation is of considerable value in 
suggesting the LDE method. Without the supersymmetric formulation, it is not obvious 
that modifying the physical Hamiltonian (p^) as in Eq. ( f2?D is a useful step. Casting a 
quantity in a different form can often suggest new methods of approximation. 

We first consider a zero- dimensional density of states model which provides support for 
the convergence of the LDE method. The model is essentially the dimensional reduction of 
expression 



dx*dxds*ds s*sexpi 



A 



(29) 



Apart from the s*s factor outside the exponential, this integral is just the supersymmetric 
version of the integral studied in Ref. O. Carrying out the integrations, we obtain 



2tt 3 /A e a2/2X erfc(a/v / 2A) 



(30) 



(Note that, without the s*s factor, the integral would equal one.) Let us now approximate 
the integral ( p9|) using the LDE method. Modifying the argument of the exponential as in 
Eq. (p6|), we have 



I s = J d X *d X ds*ds s'sexpi j-fiStS + fl (fi - a)$ f $ - ^ (V$) 2 j. 



(31) 
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The terms in the ^-expansion are simple enough that a general expression can be written 
down for the expansion up to arbitrary order N: 

We were able to evaluate In, with Q fixed by the PMS condition, up to high order in N 
for a range of complex parameter values a and A and in each case found that the sequence 



In(Qn) converged to the exact solution (|30|) 



Thus encouraged we use the LDE to approximate the d — 1 density of states. This 
example has been solved exactly using a variety of methods&ffl enabling an immediate 
check of the accuracy of the LDE. Expanding to first order in 6 either Eq. ( p4[) with S 
replaced by Sg, or Eq. ( |28| ) after disorder averaging, we obtain (with 5 set equal to 1): 

(33) 



Setting to zero the derivative with respect to Q of the term in brackets, we obtain the 
following PMS condition: 



n s/2_ EQ i/2 + ^ J|^A = 0. (34) 

This equation admits three independent solutions for Q. However, only one of these solutions 
is physical, i.e., yields a positive nonzero density of states when substituted into Eq. (j33|). 
To second order in 5, we have 



Sh 



15ft- 1 / 2 - 5MT 3 / 2 + 3E 2 n- 5 / 2 + 6m/5 AIT 2 



txTiS 2 \ V h 



-^^XE^-^n-^j. (35) 

The PMS condition now admits six independent solutions for Q of which only one is physical. 
In Fig. we plot the order 5 and 5 2 approximations to the d = 1 density of states. Also 
shown are the exact curve and the self-consistent Born approximation (SCBA) (see, e.g., 
Ref . [T3| and references therein) . A comparison between the exact curve and the order S and 
5 2 approximations provides strong evidence for the rapid convergence of the LDE method. 
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Note that even the lowest, order 5, approximation is superior to the SCBA result. The 
convergence appears to be slowest in the region of the exponentially decaying exact density 
of states tail. It is possible that a different ^-modification procedure from that used above 
can be found which gives better convergence in this region. 

The analogous calculations for d > 2 should involve little extra difficulty The only 
new feature is the divergent nature of the terms in the 5-expansion for d > 2, so that the 
procedures of regularization and renormalization are required. Other investigations^ suggest 
that these procedures are straightforward to implement within the LDE method. It should 
also be possible to use the LDE method to approximate the disorder-averaged density of 
states for random potentials satisfying other correlation relations. 



IV. THE CONDUCTANCE 



In this section we present some ideas concerning the application of the LDE method to 
the disorder-averaged conductance. Consider a channel with cross-section A and length L 
which is connected adiabatically at both ends to reservoirs. The single electron Hamiltonian 
is given by Eqs. (^) and (|23|) with an additional hard-wall confining potential restricting 
the electron to move within the cross-section A. The supersymmetric formulation of the 
zero-temperature, disorder-averaged conductance is 



dx J A dr' ± J D$\D$ 1 D$lD$s 



— e 2 h 3 r+ L / 2 r r+ L / 2 
C = o 2T2 / dx / dr± / 

6irm z L z J-L/2 J A J-L/2 

[ Sl (r)d x ,sl(r')s 2 (r')d xS ;(r) + d xSl (r)sl(r')d x ,s 2 (r')s* 2 (r) 
-9 I si(r)a x /st(r / )s 2 (r / )s2(r) - si(r)st(r / )a x /s 2 (r / )a 2; s;(r)] 



x ( e ^ + - + e ^- + + e ^ ++ + e <s~ 



where r = (x, r±) = (x,y, z), and the actions S ab are defined as follows: 

h 2 



S ab = - J dxjdrAn 



ab 
ran 



2m 



iX 



(36) 



(37) 
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with the supervector fields $ m , m — 1,2, vanishing at the cross-section boundaries. The 
field metrics are defined as follows: 

Vmn = $mni Vmn = ~^mni Vmn = ( — -0 &mm Vmn = ( — -Q ^mm Tfl^Tl = 1,2. (38) 

Eq. ([5B|) follows from the formula for the disorder-averaged conductance of a finite-length 
wire in terms of the averaged products of advanced and retarded Green functions (see, e.g., 
Ref. |20|) and the formulation of these products in terms of supersymmetric path integrals.^ 



We can decompose Eq. ( |36|) into a sum of four terms, each involving a different action. 
It is most natural to apply the LDE method to each of these four terms separately. The 
action S ++ is very similar to the density of states action fl25|) , the only difference being the 



number of supervector fields. We therefore expect that the action S ++ can be 5-modified 



with a single parameter f2 exactly as in Eq. fl26|). The series expansion in 5 is obtained 
and the PMS condition applied after all the integrals have been carried out. The complex 
conjugate of the resulting LDE approximation yields the approximation to the S term. 

The application of the LDE method to the terms involving the actions and S~ + 
is less straightforward as a consequence of the Minkowskian field metric appearing in the 
actions. As a first step, we consider a zero-dimensional conductance model analogous to 
the density of states model defined in Eq. (^9|). These models are particularly useful for 
testing the convergence of a given ^-modification. Dimensionally reducing the S~* term in 
Eq. (p6|), we arrive at the following model: 

I = J dxidxidx*2dx2ds* 1 dsids* 2 ds2 Sis\s2S* 2 exp(iS'), (39) 

where 

S = a ($j <£>! - $^$ 2 ) + ib (${ $1 + $2> 2 ) + — ($1$! - <l4<I> 2 ) 2 . (40) 

The parameters a and b are real with b > 0. We choose b to be much smaller than A. The 
term ib is the analogue of the ie term in Eq. (|37|). In order that Eq. (|39| ) not diverge, however, 
b must be positive and non-zero, rather than infinitesimal. The integrals are readily carried 
out and we obtain 
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I = ^2n 5 /X b- 1 { e (b+ia)2/2X erfc[(6 + ia)/V2X] + e (6 " ia)2/2A erfc[(6 - ia)/V2X}} . (41) 



Note that Eq. (|1]) is just twice the real part of expression fl3"0|) multiplied by the factor n/b 
and with the parameter a replaced by a + zfe. How might we 5-modify the "action" S? The 
part Sq should have some resemblance to S and so a natural choice is 



S s = ft - $£$ 2 ) + ib + $2> 2 ) 



+ 5 



(42) 



where ft = fti + zft 2 , with fix > and ft 2 > 6 > 0. Replacing S 1 by 5*5 in Eq. (|39|) and 
expanding in 5 we find, however, that the integrals are infinite since the integrands diverge 
for sis* — S2S2 — > —00. The solution to this problem is to split the range of integration in Eq. 
( P^) into the two regions Sis* > s 2 ^2 and SisJ < s 2 S2 and apply the LDE method separately 
to each integral with ft replaced by ft* = fti — ift 2 in the integral over the latter region. 
Since the two integrals are complex conjugates of each other, it is sufficient to consider only 
one of them. We are therefore led to consider the following 5-modification: 

I 5 = 2 Re J dx* 1 dxidx* 2 dx2ds* 1 dsids* 2 ds 2 Sxs^s^ 0($l$i - $ 2 $ 2 ) exp(iS s ), (43) 

where Ss is given by Eq. ([42]) and 

e($t $i _ $ t $2 ) = J_ f + °° dq ^ q _ i£ yi eX p[^($t $1 _ $t $2)] 

£1X1 J-00 

= e(sjsi - s^s 2 ) + (x*iXi - X2X2) S(slst - s* 2 s 2 ) 

-x\XiXlX2 S'{s* lSl - s* 2 s 2 ). (44) 

The function 0($J$ 1 — $2^2) limits the integration range to Sis* > s 2 S2, while preserving 
super symmetry. We call this the supersymmetric step function. Using the identity 0($|$i — 
$i>$ 2 ) + ©($5$ 2 - = 1, we can check that Is=i = I. Expanding Eq. ( f4"3|) with respect 
to 5 and then carrying out the integrals, we obtain an expression which coincides with twice 
the real part of expression (|32| ) multiplied by the factor n/b and with the parameter a 
replaced by a + ib. Thus, we can immediately conclude that applying the PMS condition 
and then taking the real part yields a series which converges to the exact solution (|4~Tf) . 



Having found a way in which to apply the LDE method to the zero-dimensional conduc- 
tance model, we can now try to apply the method in the same way to the S + ~ and S~ + 
terms in Eq. fl36|). The 5-modified actions and S$ + analogous to Eq. P^) are clear. As 
for the supersymmetric step function, the natural choice is 

9 (/ dr rC*l*n) = 2^ £7 Ml - ^ exp (iq J dv v^a) • (45) 

The delta-modifications of the terms involving S + ~ and S~ + in Eq. ( |36l) are then just the 
analogues of Eq. fl43|) . We are faced with some difficulty, however, in evaluating the terms 
in the 5-expansion series. The difficulty lies in performing the path integrals in the presence 
of the step function (|45|). The first thought is to carry out the path integrals first and then 
the g-integral appearing in the definition of the step function. But interchanging the order 
of integration in this way makes the path integrals ill-defined for the same reason as was 
mentioned above for the ordinary integrals of the zero- dimensional model. One possible 
resolution is to have a series expansion of the step function involving only polynomial and 
Gaussian functions of the superfields — so that the path integrals can be performed — and 
to link the order of this expansion with the order of the (^-expansion. Relevant discussions 
can be found in Ref. [1]. 

V. CONCLUSION 

We have applied the LDE approximation method to the quantum dynamics of a sin- 
gle electron in a random potential. The LDE method was originally developed for the 
study of quantum field systems such as A0 4 field theory. The idea that the LDE method 
might be applied to disordered systems as well comes from the resemblance of the path 
integral expressions of the A0 4 field theory n-point functions to the supersymmetric path 
integral expressions of quantities such as the averaged density of states and conductance for 
Gaussian delta-function-correlated disorder. The supersymmetric actions associated with 
these quantities contain several unusual features, such as Faddeev-Popov rather than Dirac 
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kinetic term for the anticommuting fields, imaginary coupling, and Minkowskian field co- 
ordinate metric. We showed how the LDE method can be applied to such systems using 
as illustrative examples the ground state energy of a supersymmetric anharmonic oscillator, 
a zero-dimensional density of states model, the one-dimensional density of states, and a 
zero-dimensional conductance model. 

The next stage is to apply the LDE approximation to the density of states in d > 2 
for Gaussian delta-function as well as other types of correlated disorder. The calculations 
should be similar to those for the d = 1 density of states, the only essential new feature being 
the need to regularize and renormalize. The ideas presented concerning the application of 
the LDE method to the conductance must also be developed further. 

We now finish with a brief description of another possible application of the LDE method. 
This concerns the use of a supersymmetric non-linear cr-model to determine the conductivity 
critical exponent for the Anderson metal-insulator transition (see, e.g., Sects. 3 and 4 of 



Ref. [10] and also Ref. for a review of the Anderson transition). The cr-model is usually 
quantized using the 2 + e expansion.il'i! However, the conductivity exponent for e = 1 
(d = 3) was found not to agree with the accepted value from numerical calculations.il One 
possible reason put forward for this disagreement is the omission from the cr-model action of 
high-order gradient terms which may be relevant to the fixed point and e-expansion (see, e.g., 
Ref. ^ for a review). Another possibility is that the cr-model is adequate, but the method 
of perturbative quantization is not. With respect to the latter possibility, it would be of 
interest to try to apply the LDE or some related method to the supersymmetric non-linear 
cr-model in order to determine the conductivity exponent. 
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FIGURES 

FIG. 1. Ground state energy in units of \u\ versus uj 2 /X 2 / 3 . The dashed line is the order S 
approximation and the solid line the order # 3 approximation. 

FIG. 2. Energy density of states per unit length versus Fermi energy for disorder strength 
A = 1/V2 in units h = m = 1. The solid line is the exact curve, the dashed line the order 5 approx- 
imation, the dotted line the order 5 2 approximation, and the dot-dashed line the self-consistent 
Born approximation. 
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Fig. 2 Blencowe 



exact 




